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In  this  article,  we  are  interested  in  the  performances  improvement  of  Douhly 
Fed  Induction  Machine  (DFIM)  operating  in  motor  mode  hy  the  use  of  the 
direct  torque  control  (DTC).  Firstly,  we  focused  on  the  modeling  of  the 
DFIM  and  the  study  of  the  functioning  principle  of  the  DTC  control.  Then, 
we  implement  this  control  on  the  Matlab/Simulink  environment.  Secondly, 
we  present  the  simulation  results  of  the  proposed  control.  The  analysis  of 
these  results  shows  clearly  that  the  system  based  on  the  DFIM  studied 
follows  perfectly  the  set  points,  what  allowed  us  to  justify  the  efficiency  of 
the  elaborate  control. 
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I.  INTRODUCTION 

Currently,  the  doubly  fed  induction  machine  (DEIM)  is  widely  used  in  high  power  industrial 
applications,  due  to  many  advantages  over  other  rotating  electrical  machines  [1].  The  supply  of  the  DEIM 
working  in  motor  mode  and  in  variable  speed  is  assured  by  two  inverters  of  voltage  which  are  reciprocally 
connected  to  the  stator  and  rotor  windings  [2].  The  control  of  this  machine,  which  is  essentially  non-linear 
due  to  the  coupling  between  the  flux  and  the  electromagnetic  torque,  is  relatively  complex  [3].  Several 
control  strategies  are  used  to  control  the  DEIM  to  obtain  a  decoupling  between  the  torque  and  the  flux 
[4]-[5].  But,  we  always  look  by  using  innovative  solutions  of  control  to  improve  its  performances. 

The  direct  torque  control  (DTC)  present  an  attractive  solution  to  have  a  functioning  with  better 
performances  of  this  machine,  in  the  variable  speed  applications.  This  control  strategy  was  introduced  by 
Takahashi  [6]-[7]  and  Depenbrock  [8],  it  is  mainly  characterized  by  a  good  dynamic  response  of  the  torque,  a 
good  robustness  and  a  lesser  complexity  with  regard  to  other  control. 

In  this  work,  we  present  the  direct  torque  control  of  a  DFIM  functioning  in  motor  mode.  The 
objective  is  to  guarantee  the  best  performances  regarding  the  robustness  vis-a-vis  the  variations  of  its 
parameters  and  the  load  torque.  The  first  part  of  this  article  will  be  devoted  to  the  modeling  of  the  DFIM.  The 
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second  part  deals  with  the  control  of  this  machine.  In  a  third  part,  in  order  to  validate  our  model,  simulation 
results  using  the  Matlab/Simulink  software  will  be  presented  and  analyzed. 


2.  MODELING  OF  THE  DFIM 

The  global  system  is  illustrated  in  Figure  1  below,  it  consists  of  two  converters,  one  is  connected  to 
the  stator  and  the  other  is  connected  to  the  rotor  of  the  machine  studied  [9] -[10]. 


Stator 


Figure  1.  Overall  scheme  of  the  system  studied 


The  model  of  the  doubly  fed  induction  machine  after  the  Park  transformation  is  defined  by  the 
electrical,  magnetic  and  mechanical  following  equations  [11]-[12]: 

2.1.  The  electrical  equations 

The  electrical  stator  and  rotor  equations  of  the  DFIM  in  the  reference  frame  (d,  q)  are  expressed  by 
[13]-[14]; 


Vsi-Rs-isi+^-^-Vsc, 

^sq  -  Rs-  isq  +  ■  VsH 

V  —Pi  -I-  m 
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2.2.  The  magnetic  equations 

The  expressions  of  the  flux  are  extracted  from  the  electrical  Equations  (2): 


(Vsd  ~  ^S‘  ^sd  F  M.  tj-(j 
Vsq  ~  ^S‘  ^sq  F  M.  i-irq 

I  *Prd  ~  ^r-  Wd  F  M. 

<Vrq  ~  ^r-  Wq  F  M.  igq 


(3) 


2.3.  The  mechanical  equation 

The  fundamental  equation  of  dynamics; 


Tem  =  Tr+J.^  +  f.il 


(4) 
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The  electromagnetic  torque  is  expressed  as  a  function  of  the  currents  and  the  flux  by: 

Tem  ~  P-  {sPsd-  ^sq  T  Vsq-  ^sd)  (5) 


With: 

Vs(d,q),  Vr(d,q) 

:  Stator  and  rotor  voltages  in  the  reference  frame  (d,  q). 

Is{d,q],  Ir(d,q) 

:  Stator  and  rotor  currents  in  the  reference  frame  (d,  q). 

tPs(d.q),  tpr(d.q) 

:  Stator  and  rotor  flux  in  the  reference  frame  (d,  q). 

Rs,  Rr 

:  Stator  and  rotor  resistances. 

Est  Er 

:  Stator  and  rotor  inductances. 

M 

:  Mutual  inductance. 

P 

:  Pole  pair  number. 

a 

:  Dispersion  coefficient. 

;  Stator  and  rotor  pulsations. 

Tern,  Tr 

:  Flectromagnetic  and  load  torque. 

Q 

:  Rotation  speed  of  the  machine  (a)=p.Q). 

J 

;  Moment  of  inertia. 

fi 

:  Coefficient  of  viscous  friction. 

The  implementation  on  Matlab/Simulink  of  the  DFIM  is  represented  by  the  following  Figure  2; 


Figure  2.  Simulink  Model  of  the  DFIM 


3.  DIRECT  TORQUE  CONTROL 
3.1.  Principle  of  the  DTC  control 

The  general  principle  of  direct  torque  control  is  based  on  the  direct  regulation  of  the  torque  and  flux 
of  the  doubly  fed  induction  machine  by  selecting  one  of  the  eight  voltage  vectors  generated  by  each  voltage 
inverter.  This  choice  is  generally  based  on  the  use  of  hysteresis  regulators  whose  function  is  to  control  the 
state  of  the  machine  [15]- [16].  The  voltage  vector  V  can  be  written  in  the  form  [17]; 

F  i2ii  i^n 

F  =  l/o(5a  +  S^eT-  +  S,e^)  (6) 

o  Si(i=a,b,c)=0  if  the  high  switch  is  closed  and  the  low  switch  is  open, 
o  Si(i=a,b,c)=l  if  the  high  switch  is  open  and  the  low  switch  is  closed. 

The  combinations  of  three  magnitudes  (Sa,  Sb,  Sc)  allow  to  generate  eight  positions  of  the  voltage 
vector  V  which  two  correspond  to  the  zero  vector  (Sa,  Sb,  Sc)  =  (000)  and  (1 1 1)  as  shown  in  Figure  3. 
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Figure  3.  Elaboration  of  voltage  vectors  from  voltage  inverters  of  the  DFIM 


3.2.  Strategy  of  the  DTC  control 

To  realize  this  command  it  is  necessary  every  instant  to  control  the  stator  flux,  the  rotor  flux  and  the 
electromagnetic  torque  of  the  machine. 

3.2.1.  Control  of  the  stator  and  rotor  flux 

The  relation  between  the  flux  vectors  of  the  DFIM  and  the  voltage  vectors  is  given  by  the  following 
equations  system  [18]: 

1^(0  =  /o  (l?  -Rs-^-dt  +  $so 

=  /o(i?  -  Rr  X)- dt  +  $ro 


On  an  interval  [0,  te],  we  consider  that  the  terms  RJs  and  RJr  are  respectively  negligible  in  front  of 
Vs  and  Vr.  Thus,  we  can  write  the  equations  (7)  in  the  form: 

f  0s  ~  0sO  F  Vs-te 
0rO  F  Vf.  tg 

In  general,  for  each  flux  we  have: 

0(k+l)  =  0(fc)  +  ^-ie  (9) 


With: 

;  Vector  of  the  flux  in  the  step  of  current  sampling. 

0(fc+i)  ■  Vector  of  the  flux  in  the  step  of  following  sampling. 

tg  :  Sampling  period. 

In  order  to  obtain  a  functioning  with  a  module  of  constant  flux,  we  can  follow  the  extremity  of  (p  on 
an  almost  circular  trajectory,  by  choosing  a  correct  sequence  of  the  vector  V  for  the  different  successive  time 
intervals  of  duration  4. 


Figure  4.  Example  of  the  evolution  of  (p  extremity 
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3.2.2.  Control  of  the  electromagnetic  torque 

The  expression  of  the  electromagnetic  torque  in  the  reference  frame  (a,  P)  is  written  [19]: 

=  (10) 

With:  S  =  Qg  —  By-,  the  angle  of  phase  shift  between  the  two  stator  and  rotor  flux. 

According  to  the  equation  (10),  we  notice  that  the  electromagnetic  torque  of  the  DFIM  depends  on 
the  amplitude  of  two  vectors  (p^  and  (p-^  and  the  angle  between  them,  if  we  succeed  in  controlling  perfectly 
this  two  flows  in  module  and  in  position,  we  can  thus  control  the  torque. 

3.2.3.  Choice  of  the  voltage  vector 

The  voltage  vector  selection  depends  on  the  flux  position  in  the  reference  frame  (a,  P)  and  the 
variation  desired  for  the  flux  module  and  for  the  electromagnetic  torque.  The  evolution  space  of  the  flux  is 
decomposed  into  six  sectors  of  60  °:  When  the  flux  vector  is  in  the  sector  Si,  the  control  of  the  flux  and  the 
torque  can  be  made  in  the  way  presented  in  the  following  Table  1: 


Table  1.  Choice  of  the  voltage  vector 


The  evolution  of  the  flux  and  the  torque 

Choice  of  the 

Flux 

Torque 

voltage  vector 

Increasing 

Increasing 

Vi+i 

Decreasing 

Increasing 

Vi+2 

Increasing 

Decreasing 

V,, 

Decreasing 

Decreasing 

Vi.2 

The  vectors  Vo  and  V?  are  chosen  to  stop  the  movement  of  the  flux. 

3.3.  Estimation  of  the  flux  and  the  torque 

The  estimation  of  the  stator  and  rotor  flux  is  basing  of  the  stator  and  rotor  voltages  and  currents.  The 
flux  components  of  the  DFIM  in  the  reference  frame  (a,  P)  are  determined  by  [20]: 

^sa  ~  /q  (Vsa  ~  ■  ^sa')' 

Vsp~  /q  (^s/?  ~  4s^).  dt  (11) 

^ra  ~  Jq  C^ra  ~  ■  Wa)- 

SPrp  ~  /q  (W/?  “  B.J-  .  irp')-  dt 


With:(/)^  =  +  j-  <Ps/?  4’r  ~  Vra  +  j-  The  modules  of  the  stator  and  rotor  flux  are  written: 


0s  ^sa  T  ^s/? 

I - 

0r  ~  J^ra  T 


The  positions  of  the  flux  vectors  are  calculated  as  follows: 


f0,^arctg(^) 
=  arctg  (^) 

^  H'ra 


(13) 


The  electromagnetic  torque  is  estimated  from  the  measured  stator  currents: 

Tern  ~  P'  (*Psa'  ^s/?  ~  *Ps/?'  ^sa)  (14) 

The  DTC  control  is  vectorial,  it  is  necessary  to  decompose  three  currents  and  voltages  (iabc,  Vabc)  in 
components  direct  (ia,  Va)  and  quadratic  (ip,  Vp)  according  to  the  Concordia  transformation,  such  as: 
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(Xa 

\Xp, 

i-Ji) 

1  -i 

2 

0  ^ 

3 

1  1 

(1) 

(15) 

With:  X  = 

Vr,  is  or  ir 

3.4.  Elaboration  of  the  control  vector  (Switching  table) 

3.4.1.  Regulators  elaboration  of  the  flux  and  the  torque 

The  difference  e^,  between  the  reference  flux  <j)ref  and  the  estimated  flux  <j)est  is  introduced  into  a 
hysteresis  regulator  with  two-level,  the  latter  will  generate  the  binary  variable  {H^  =  0,  1)  at  its  output, 
representing  the  desired  evolution  for  the  flux: 

a.  7/0=0  means  that  the  flux  must  be  decreased. 

b.  7/0=1  means  that  the  flow  must  be  increased. 


Figure  5.  Hysteresis  regulator  with  two  levels 


Similarly,  the  immediate  error  of  the  couple  STem  is  calculated  by  comparing  its  reference  value  Tem-ref 
and  its  estimated  value  Tem-est,  then  applied  to  a  three-level  hysteresis  regulator,  generating  at  its  output  the 
variable  77rem  with  three  levels  (-1,0,1),  representing  the  direction  of  temporal  evolution  desired  for  the  torque: 

a.  77rem=  1  means  that  the  torque  must  be  increased. 

b.  77rem=  0  means  that  the  torque  must  be  decreased. 

c.  77rem=-l  means  that  it  is  necessary  to  maintain  the  torque  inside  the  band. 


Figure  6.  Hysteresis  regulator  with  three  levels 


3.4.2.  Elaboration  of  the  switching  table. 

The  state  choice  of  each  voltage  inverter  is  made  in  a  switching  table  elaborated  according  to  the 
state  of  the  variables  (77^  and  Hrem)  at  the  output  of  the  flux  regulators  and  the  electromagnetic  torque,  as 
well  as  sector  giving  the  information  about  the  position  of  vector  of  the  flux.  Table  2  is  presented  in  the 
following  form: 


Table  2.  Switching  table 


Hi,  ox  Hi, 

tiTem 

1 

2 

Sector  Si 

3  4 

voltage  vector 

5 

6 

1 

1 

V2(110) 

VjCOIO) 

V4(011) 

VsCOOl) 

VsClOl) 

V,(100) 

1 

0 

V7(lll) 

v„(000) 

V7(lll) 

Vo(OOO) 

V7(lll) 

v„(000) 

1 

-1 

V6(101) 

V,(100) 

V2(110) 

VjtOlO) 

V4(011) 

VsCOOl) 

0 

1 

VstOlO) 

V4(011) 

VstOOl) 

V6(101) 

Vi(lOO) 

V2(110) 

0 

0 

V„(000) 

V7(lll) 

Vo(OOO) 

V7(lll) 

Vo(OOO) 

V7(lll) 

0 

-1 

VsCOOl) 

Ve(lOl) 

Vi(lOO) 

V2(110) 

VjCOIO) 

V4(011) 
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The  general  structure  of  direct  torque  control  of  the  doubly  fed  induction  machine  is  detailed  in  the 
following  Figure  7 : 


Figure  7.  General  structure  of  the  direct  torque  control  applied  to  the  DFIM 


4.  RESULTS  AND  DISCUSSION 

The  global  model  of  the  system  studied  (DFIM  and  blocks  of  the  DTC  control)  is  simulated  in  the 
Matlab/Simulink  environment. 
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— 

— 
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Figure  8.  Simulation  scheme  of  the  DTC  control  on  the  Matlab/Simulink  environment 


4.1.  Simulation  results 

To  show  the  performances  of  direct  torque  control  applied  to  DFIM,  we  present  the  simulation 
results  of  the  system  for  reference  flux  (cps-ref  =  IWb,  cpr-ref  =  0.5Wb).  The  bandwidth  of  hysteresis 
comparator  of  the  torque  is  fixed  in  +0.01N.m  and  those  of  the  comparators  of  the  stator  and  rotor  flows  are 
fixed  to  +0.005  Wb. 
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4.1.1.  Test  with  constant  speed 

The  figures  below  present  the  response  of  the  magnitudes  of  the  doubly  fed  induction  machine  at 
constant  speed  and  variable  load  (the  reference  torque  of  set-points  ranging  from  0  — >  10  N.m  — >  5  N.m).  In 
this  test,  the  rotation  speed  reaches  its  reference  without  overtaking,  it  is  insensitive  to  changes  in  the  load 
(Figure  9. a).  Figure  9.b  show  that  the  torque  has  good  dynamics  whose  mean  value  follows  the  setpoint 
values  with  a  very  fast  response  time  which  equals  t,  =  0.32s,  but  it  has  low  ripples.  The  stator  and  rotor  flux 
remain  constant  at  its  reference  values  whatever  the  applied  load  which  shows  that  the  torque  and  the  flux  are 
decoupled  (Figure  9.c).  The  evolution  of  these  flows  in  the  reference  frame  (a,  P)  is  perfectly  circular 
(Figure  9.d).  The  components  of  the  stator  and  rotor  currents  in  the  (a,  b,  c)  reference  have  sinusoidal  patterns 
with  a  frequency  proportional  to  the  reference  speed,  they  reply  well  to  the  variations  imposed  by  the  load 
torque  (Figure  9.e  and  Figure  9.f). 


1,5 

Time  (s) 

(a) 


(c) 


Figure  9.  Test  with  constant  speed,  (a)  Rotation  speed  (b)  Electromagnetic  torque  (c)  Stator  and  rotor  flux  (d) 
Evolution  of  the  stator  and  rotor  flux  in  the  reference  frame  (a,  P)  (e)  Stator  currents  (f)  Rotor  currents 


4.1.2.  Test  with  variable  speed 

The  figures  below  show  the  simulation  results  of  the  DTC  control  applied  to  DEIM  with  variable 
speed  and  variable  load.  In  the  variable  speed  test,  the  Eigure  10. a  and  Eigure  lO.b  show  the  high  dynamics 
of  the  rotation  speed  and  the  electromagnetic  torque  of  the  DEIM,  with  a  fast  response  time  to  the  change  of 
the  set-points.  The  components  of  the  stator  and  rotor  currents  in  the  (a,  b,  c)  reference  are  always  sinusoids, 
their  frequencies  and  their  amplitudes  vary  respectively  as  a  function  of  the  rotation  speed  and  the  load 
torque  imposed  (Eigure  lO.e  and  Eigure  lO.f).  The  stator  and  rotor  flux  remain  constant  at  its  reference 
values,  which  validates  that  the  DTC  control  is  robust  towards  the  variation  of  the  load  and  the  speed  as 
shown  in  Eigure  lO.c.  These  stator  and  rotor  flux  are  trapped  in  a  circle  regardless  of  the  applied  load,  which 
shows  that  the  flux  constancy  despite  the  load  as  shown  in  Figure  10. d,  we  conclude  that  the  torque  and  flux 
are  decoupled.  These  simulation  results  show  the  effectiveness  of  this  control  technique. 
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Figure  10.  Test  with  variable  speed,  (a)  Rotation  speed  (b)  Electromagnetic  torque  (c)  Stator  and  rotor  flux 
(d)  Evolution  of  the  stator  and  rotor  flux  in  the  reference  frame  (a,  P)  (e)  Stator  currents  (f)  Rotor  currents 


5.  CONCLUSION 

In  this  work,  we  presented  the  direct  torque  control  applied  to  a  DFIM  used  in  motor  mode, 
powered  by  two  converters,  one  on  the  stator  and  the  other  on  the  rotor.  After  a  study  of  the  modeling  of  the 
system,  we  realized  the  implementation  of  the  DTC  control  based  on  the  direct  regulation  of  flux  and  torque. 
The  interpretation  of  the  obtained  results  shows  clearly  that  the  system  follows  perfectly  the  reference 
values  and,  by  consequence,  supplies  a  good  static  and  dynamic  performance  of  the  studied  machine. 


APPENDIX 


Table  3.  Parameters  of  the  DFIM 


Variable 

Symbol 

Value  (unit) 

Nominal  power 

Pm 

1.5  kW 

Stator  nominal  voltage 

V„ 

220V 

Rototor  nominal  voltage 

Vm 

BOV 

Frequency 

f 

50  Hz 

Pair  pole  number 

P 

2 

Stator  self  inductance 

Ls 

0.295  H 

Rotor  self  inductance 

Lr 

0.104  H 

Maximum  of  mutual  inductance 

M 

0.165  H 

Stator  resistance 

Rs 

1.75  £2 

Rotor  resistance 

Rr 

1.68  Q. 

Total  viscous  frictions 

f 

0.0027  Kg.mVs 

Total  inertia 

J 

0.01  Kg.m^ 
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